Is there a true Model-D critical dynamics? 
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We show that non-locality in the conservation of both the order parameter and a noncritical 
density (model D dynamics) leads to new fixed points for critical dynamics. Depending upon the 
parameters characterizing the non-locality in the two fields, we find four regions: (i) model-A like 
where both the conservations are irrelevant (ii) model B-like with the conservation in the order 
parameter field relevant and the conservation in the coupling field irrelevant (iii) model C like where 
the conservation in the order parameter field is irrelevant but the conservation in the coupling field 
is relevant, and (iv) model D-like where both the conservations are relevant. While the first three 
behaviours are already known in dynamical critical phenomena, the last one is a novel phenomena 
due entirely to the non-locality in the two fields. 



Based on the long time behaviour, the dynamics of dif- 
ferent critical systems are classified in various dynamic 
universality classes. These classes are characterized by 
several factors like the number of conserved field, non- 
vanishing Poisson brackets etc, in addition to the usual 
quantities necessary for equilibrium. Among the models 
which have no propagating modes, Model A is a system in 
which there are no conserved fields, Model B is a system 
with a conserved order parameter while Model C is the 
class of systems with a non-conserved order parameter 
field coupled to a conservative field. A fourth classifica- 
tion was also made - Model D, in which both the order 
parameter field and the coupled field are conserved. This 
model, however, was shown to be equivalent to Model B 
as the coupling to the field becomes irrelevant in the large 
length long time scale limit. 

The dynamical behaviour in systems with one or more 
conserved field like Models B, C and D in the above clas- 
sification is usually studied by considering local conser- 
vation. However, conservation is by no means necessarily 
local and in fact, non-local conservation is a more gen- 
eral condition. A few recent studies on critical or near- 
critical dynamics show that non-locality in the conserved 
field may change the dynamical class, besides the obvious 
modifications in nonuniversal details. This kind of effect 
of non-locality in the conservation has been noted in the 
study of critical dynamics of model B and C and in the 
numerical study of the early time effect of a model be- 
lieved to belong to the class of model C and on the zero 
temperature dynamics (phase-ordering kinetics or Ost- 
wald ripening ) §, |, § [||, § . Besides the theoretical 
interest on its own right, one of the important applica- 
tions of nonlocal dynamics is in speeding up numerical 
simulations, prompting implementation of various types 
of nonlocal moves Hence the need of a general classi- 
fication of effects of nonlocality in dynamics. 

In a sense, model D is the most generalised class of 
dynamical critical phenomena with dissipative dynam- 
ics. However, with local conservation, it is equivalent to 



Model B showing no novel features. Motivated by the 
results of model C j|, we investigate the problem of non- 
local conservation in Model D in the most general way 
by keeping both the conservations in the order param- 
eter field and the coupled field non-local. One of our 
main results is to show the emergence of a true model D 
behaviour by nonlocal dynamics, unlike the local conser- 
vation case. 

We introduce non-locality by using a non-local kernel 
|^|, ||, [5| in the chemical potential keeping the continuity 
equation intact. Using cj> for the order parameter and m 
for the secondary coupling parameter, the Hamiltonian 
of model D is given by 

H = J d d x[ir0 2 (x) + i(V0(x)) 2 + £0 4 (x) 

+70 2 (x)m(x) + icrWfx)], (1) 

where 2 (x) = £?=i #(x), 0*(x) = E£=i#(*)]V = 
is the mean field critical point, and C > ensures 
noncriticality of m. There is a short distance cutoff which 
in momentum space is an ultraviolet cutoff A. We shall 
set A = 1 in the calculation. The above Hamiltonian, 
Eq. (0) is the same as in Model C[Q, |). Making 7 = 
effectively reduces the above model to model A. 

The dynamics ^t] [ s given by the following equations: 
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(2b) 



[ft] Additional torques for n = 3 (Heisenberg model) may lead to 
new dynamic universality classes H. We do not consider such 
effects. 
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The noises 77 and £ obey 



(77) = 0, (C) = 0, (3a) 
Mk, *Mk', 0> - 2(-i- + i)5(k + k')*(t - 0, (3b) 
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— )<J(k + k')<S(i-i')-(3c) 



Here, k = k |, and we have introduced two different 
parameters p and a to denote the non-localities in the 
order parameter field and the coupling field respectively. 
The dissipative or the kinetic coefficients (ro and Ao) 
are denoted by the subscript while the transport coef- 
ficients (T p and A CT ) by the corresponding powers of k. 
For the local case, both p and a assume a value equal 
to 2. Results are known for some limiting values of p 
and a: a) (p = 0,7 = 0) corresponds to model A [0, b) 
(p = 2, 7 = 0) and (p = 2, a = 2) give model B like be- 
haviour L c) any value of p > and 7 = is a non-local 
Model B 2 , H , and d) p = with any c > is non-local 
model C §. 

A Hamiltonian of the type of Eq. [I] occurs near the 
critical wings of a tricritical system or in critical systems 
with additional constraints. Model D dynamics are rele- 
vant in such cases, for example in three or four compo- 
nent mixtures or in spin-1 or Potts models Jnj. A spin-^ 
example, though a bit unrealistic, would be an antifer- 
romagnetic system with conserved magnetisation (as in 
model C [|[ [ll] ) as well as conserved staggered magneti- 
sation. The Kawasaki dynamics can be used to conserve 
the magnetisation maintaining at the same time that the 
exchanges between opposite spins have to be done in such 
a way that the order parameter is unchanged. 

A momentum-shell renormalization group approach is 
used to study the long distance long time behaviour. In 
this approach, small scale fluctuations in space (between 
Ae~ si and A in momentum space) are integrated out, 
and the cutoff rescaled to A. One then obtains an ef- 
fective Hamiltonian and effective equations for the dy- 
namics valid for longer scales. The universal features are 
obtained from the fixed points of the flow equations for 
the various parameters of the problem. As per the usual 
RG scheme, any parameter that vanishes (grows) in the 
long scale limit is called an irrelevant (relevant) variable. 
A fixed point is then characterized by the set of relevant 
parameters. In this scheme the four dynamics classes 
correspond to the relevance and irrelevance of A^ and T p 
as mentioned in the abstract. 

The statics of this system is the same as in Model C 
and the fixed point values of the parameters r, u and 7 2 C 
are as given in ref B. The dynamical equations, on the 



other hand, take the forms 
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(4a) 
(4b) 
(4c) 
(4d) 



where Q = V 1 



„ , To" 1 , P = A" 1 + A" 1 and K d = 
2 1 -<* 7r -<V 2 r(d/2). The parameters A CT , A , T p , T , 7 and 
C appearing in the above equations are all I dependent. 
To simplify notation, this Z-dependence has been sup- 
pressed. 



From (4b) we get 

r (1 + x) ni 2 K d 
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where x = -E 2 - is the dimensionless ratio of the two coef- 
ficients for the primary order-parameter (note: A = 1). 
This quantity x satisfies 
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(p — Y)x, where Y 
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Here we have introduced the ratios of the transport and 
kinetic coefficients of m with the kinetic coefficient of <fi, 
namely, 
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and p = p a + po- 



By using the fixed point value of C and Eq. 
be rewritten as 



po can 



Mo 



(x + \)a/v 
2{z - ajv)' 

The flow equation for p a is 

?tl =IIa { (r -2 + a/i;-Y}. 



(7) 



We consider the different solutions of the equations for 
x and pa- Corresponding to the relevance of the con- 
served quantities, one gets different regions in the p — a 
plane with different dynamic critical behaviour. It should 
be mentioned here that for n > 4, the coupling 7 scales 
to zero and one effectively gets a non-local Model B. This 
has already been considered and shown to be model A- 
like when globally conserved ||. Our discussions are for 
the n < 4 region where the coupling survives. 

There are three fixed points of x: x = 0, x = 00 and a 
non-zero finite fixed point value for x from Eq. (^|): 
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Case 1. x — in the long length scale limit: 

This corresponds to non-local Model C: the conserva- 
tion of the order parameter is irrelevant. The fixed points 
and their stabilities of non-local model C, which were ob- 
tained in Ref . Q from the solutions of p G , now depends 
on the value of p. In general, x — will be a stable fixed 
point as long as p < Y. We focus on the stability of the 
various regions for non-zero values of p. 

(a) p a has a finite fixed point value and z — o + a/v. 
This is a stable solution in the region (2/n) — 1 > p > — 1 
(where p = for p = 0. For non-zero p this remains 
stable for p < p + 1 . 

(b) p a — oo : this is valid for p > (2/n) — 1. Here 
2 = 2 + 2a/nv. This solution a; = 0, p a — oo is stable as 
long as p < 2a /nv. 

(c) For Pa —— . the stability of x — is valid only for 
p < — 1 and p < 0(e 2 ). Here z = 2 + 0(e 2 ) and the 
behaviour is model A like. Any non-zero value of p thus 
destroys this model A like region. Therefore, a region 
where both the conservations are irrelevant is restricted 
strictly to p = 0. 

Case 2. A finite non-zero value of x in the long length 
scale limit: This occurs when p — Y. In order that Tq 
reaches a finite fixed point value, we must have z = 2 + p. 
In this case, p a again has the following fixed points: 

(a) p a = oo for a — 2 + a/v > p 

(b) pa = for a — 2 + a/v < p 

(c) pa nonzero finite for a — 2 + a/v = p 

Physically a finite valued fixed point of x means that 
the conservation in the order parameter <f> is relevant. If 
at the same time p a attains a non-zero fixed point value, 
the conservation of the coupling field m is also relevant. 
This occurs if the fixed points (a) and (c) are stable. 

For case (a), x = P~ 2a J vn anc [ f or non-negative value 
of x, p > 2a /nv. Hence in between p — 2a /nv and 
p = a — 2 + a/v, a stable region is obtained where both 
the conservations are relevant. 

For case (b), x + 1 = ^t^t^tt-. This has a finite 
valued solution for x only for p ~ 0(e 2 ). Hence for p > 
a — 2 + a/v, x — ► oo is the only solution for p = 0. It is 
a conventional model B fixed point. 

Case (c): Here z = a + a/v and in principle along the 
entire line a — 2+a/v = p a solution exists where both the 
conservations are relevant. However, there is a restriction 
on the value of a from the condition of existence of a non- 
zero finite value of p a ■ We find that 



Pa + PO 



a/v)(x + 1) 



a/v + z*-(x +1) 



(8) 



Hence with a = 2 + pa/v; we get the condition (1 + 
p) < 2(x + l)/n. Since x is also non-zero, this implies 
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FIG. 1: The different regions in the non-local Model D with 
n = 1 is shown in the rescaled p — p plane where p = p/(a/v) 
and p — „7y ■ Regions I and II are model C like with z = 
a + a/v and z = 2 + 2a/ v respectively. Region IIP is model 
A like (as in model C of Ref. 4) with z = 2 + 0(e 2 ). Reg IV 
is the "Model D" region, z = p + 2; Reg V is Model B like, 
z = p + 2. The line p+1 = p is also model D-like upto point S 
(shown by the dashed line) . There are no discontinuities of z 
along any boundary except at the boundary between IIP and 
V. 



In the rest of the p — a plane, the only solution is 
x = oo and 1/A CT = for which only the conservation of 
the order parameter is relevant. Here z — p + 2. 

All the above possibilities of different dynamic be- 
haviours are summarized in Fig. 1 in the p — a plane (for 
n = 1). The possibility of seeing a region where both 
the conservation conditions are relevant raises new issues 
like early time effect and boundary or surface effects Q 
in the new regime. These remain to be studied. 

To summarize, Model D can be viewed as the most 
general dynamic model with purely dissipative dynamics 
and we indeed obtain all the four types of behaviours as 
soon as non-locality in the conservations (as in Eq. (||)) 
are introduced. Since our general model subsumes the 
previously studied nonlocal model C and local model D, 
it is natural to expect A-like, B-like and C-like regimes. 
The model A like region, however, turns out to be unsta- 
ble for any finite value of p. The model C like regions are 
stable for small nonzero values of p but ultimately disap- 
pear for larger values of p. The most non-trivial result, 
for general values of p and a, is the appearance (Region 
IV in Fig. 1) of a region with a new dynamical behaviour 
where both the conservations are relevant. This we call 
a true model D like region - a consequence of nonlocal 
conservation laws. 
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that this solution with both x and p a finite is valid for 
p + 1 = p/(a/v) with at least p < (2/n) — 1. 
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